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S  U  i!  II  ARY 


A  method,  is  presented  whereby  the  'Slender  Body 
Theory'  can  be  applied  to  the  determination  of  the  unsteady 
aerodynamic  forces  acting  on  slender  wings  and  wing-body 
combinations  experiencing  harmonic  deformations  in  a  com¬ 
pressible  flow.  The  analysis  holds  for  subsonic 
and  supersonic  speeds,  subject  to  restrictions  which  are 
stated  and  discussed* 

A  simplification  of  the  method  is  also  introduced 
which  is  applicable  to  many  practical  cases  and  calculations 
are  performed  on  this  basis  which  lead  to  numerical  results 
for: 

1.  'Equivalent  Constant  Derivatives'  for  a  deforming 
slender  delta  wing  using  modal  functions  which  are 
polynomials  of  the  spanwise  parameter: 

2.  '  'Rigid'  Force  Coefficients  for  a  pitching  and 

plunging,  slender,  wing-body  combination. 

These  results  are  given  as  closed  expressions  and 
in  tabular  form  and  some  of  the  results  are  also  shown  in 
graphical  form. 

Both  the  derivatives  and  the  'Rigid'  force  coeff¬ 
icients  are  defined  in  such  a  way  as  to  agree  with  the  usual 
British  Sign  Convention. 
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LES'J  OF  SYMBOLS 


The  use  of  a  ’bar'  over  a  symbol  denotes  that  it 
is  the  amplitude  of  the  harmonically  oscillating  quantity 
represented  by  the  symbol  itself. 

Symbol  Description  Defined  in 

”  ™  Section 


A,B 

P(x) 

Gs 

H(8) 

L 

V  LS 

V  L& 

M 

V  Ma 

11  (X) 

Ne^2^ 

n 

0(  ) 

Pn 

«r 

R(x) 

Ro 

S 

T 


portmanteau  symbols  in  eqn.  to  Appendix  I 

reference  axis 

function  of  II (x)  7 

constants  in  H(8)  4 

torsional  mode  4 

unsteady  lift  on  rigid  wing  or  body  6 

’Rigid’  force  coefficients  6 

constants  in  series  expansion  for  0  3 

unsteady  moment  on  rigid  wing  or  body  6 

'Rigid'  force  coefficients  6 


freestream  Mach  number  2 

modified  Mathieu  function  of  the  third 
kind  3 

'of  order' 

constants  in  series  expansion  for  0  3 

Lagrangian  generalised  force  4 

local  radius  of  body  of  revolution  7 

Maximum  '  '  '  '  '  7 


planform  area  of  wing  Appendix  I 

factor  giving  length  of  wing-body 
combination  forebody  as  proportion 
of  wing  root  chord  Appendix  VI 
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u,v,w 
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x,y,z 

zo 

r  (x) 

A 


ii 

a 


/9 

Trs,trs,Grs 

6 

e 

X> 

r\ 

K 

x 

v 


P>  P, 


CK? 


b/2 

Mathieu  function  (periodic) 
time 

perturbation  velocities 

with  suffices  -  various  upwash 
conditions 

right-handed  Cartesian  coordinates 

equation  of  reference  axis 

amplitude  at  reference  section  in 
flexural  mode 

local  cross-sectional  area  of  body- 

factor  giving  length  of  conical  nose 
on  body- 


amplitude  at  reference  section  in 
torsional  mode 

functions  of  X 

flutter  force  coefficients 

non-dimensional  spanwise  parameter 

dimensionless  amplitude  or  thickness 

elliptical  coordinate J  polar  angle 

elliptical  coordinate 

_  v  _  & 

a  “  s 

delta  planform  factor 

angular  frequency- 

local  and  freestream  densities 
R 

=  ~  ratio  of  max.  body  radius  to 
o  max.  wing  span 


3 
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7 
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perturbation  velocity  potential 


co.co  .03  local,  root  and  mean  frequency 

'  r'  m 

parameters 


w 


VC 


etc.  Suffix  c  is  used  for 
wing-body  combination 
in  place  of  *r 1 
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Appendix  V 


1 .  Introduction 


In  this  paper  a  method  is  given  whereby  the  aero¬ 
dynamic  forces  can  be  calculated  for  slender,  low  aspect 
ratio  wings  deforming  harmonically  in  a  compressible  flow. 

The  method  is  applied  to  a  slender,  cropped,  delta 
wing  and  certain  flutter  modes  are  assumed  which  take  the 
form  of  polynomials  in  the  spanvri.se  parameter.  Freedom 
of  the  wing  root  is  allowed  for  so  that  body  freedoms  can 
be  included.  In  the  latter  part  of  the  paper  the  aerodyn¬ 
amic  forces  on  a  pitching  and  plunging,  slender,  wing-body 
combination  are  evaluated. 

The  basis  of  the  method  is  the  * Slender-Body  Theory' 
which  has  been  applied  in  connection  with  the  (quasi-steady) 
stability  derivatives  for  slender,  wings  and  wing-body  com¬ 
binations  (refs.  '1,2,3,4j5>6). 

The  application  to  an  oscillating  and  deforming 
wing  has,  very  recently,  been  studied  by  Merbt  and  Landahl 
(ref,  8). 


The  solution  of  the  'cross-sectional’  problem, 
for  the  wing,  is  analogous  to  that  of  a  two-dimensional 
flat  plate  oscillating  in  a  compressible  flow  and  has  been 
treated  by  Tirnman*  (ref.  9)  and  Reissner  (ref.  10) . 

The  use  of  the  *  Slender-Body  Theory’  allows  the 


*  Only  the  regular  part  of  their  solution  is  required 
in  this  case. 


analysis  to  apply  at  subsonic,  -  and  supersonic 

speeds  subject  to  certain  restrictions  on  Aspect  Ratio,  Mach 
number,  Frequency  Parameter,  Slenderness,  Validity  at  subsonic 
Speeds  depends  on  an  approximate  satisfaction  of  the  Kutta- 
Jowkowski  condition. 

The  assumptions  of  linearised,  thin  aerofoil  theory 
are  used,  the  fluid  is  perfect,  the  flow  irrotational  and 
harmonic  motions  are  considered  throughout. 


2,  The  Slender-Bodv  Theory 

The  coordinate  system  used  is  shown  in  Figure  1 
where  right-handed  rectangular  axes  are  drawn  from  an  origin, 
0,  fixed  in  the  wing,  with  the  x-axis  parallel  to  the  main 
stream  and  the  z-axis  upward.  It  is  assumed  that  the  wing 
is  a  thin,  flat  plate  oscillating  about  its  position  of  zero 
incidence  in  the  plane  z  =  0,  but  always  lying  in  the 
immediate  vicinity  of  the  plane. 

The  perturbation  velocity  potential,  0,  satisfies 
the  equation 

v2*  =  j  (vfe  ♦  I if  *  . ^ 

The  conditions  holding  at  the  surface  of  the  wing 
are  specified  by  a  prescribed  'downwash'  w(x,y,0,t)  and 
the  stipulation  that  the  relative  normal  velocity  of  the  air 
and  of  the  wing  is  zero, 

Applying  the  assumptions  of  the  'Slender  Body 
Theory'  implies  that  the  x-derivatives in  equation  (2,1  J  are 
neglected  and  the  two-dimensional  flow  at  any  cross-section 


is  then  given  by  the  wave  equation, 

^yy  +  ^zz  =  T  *  ^tt  . (2*2) 

a 

The  approximation  equation  (2,2)  is  satisfied 
(ref.  1J,5)  if: 

11 -4J  (£f«'  . ^ 

where  is  the  freestream  Mach  number  and  s  is  the 
local  semi-span,  or  for  a  triangular  wing  if 

|l  “  *4,|  AR2  <<16  . . (2.3)1) 

where  JR  is  the  aspect  ratio. 
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If  the  influence  of  the  time-derivative  tenn  in. 
(2.2)  is  so  small  that  it  can  be  neglected  the  equation 
reduces  to  Laplace* s  equation,  as  for  steady  flow,  and  the 
difference  between  the  unsteady  and  steady  flow  cases  mani¬ 
fests  itself  entirely  in  the  linearised  Bernoulli  Pressure 
equation, 


This  implies  that  the  root  frequency  parameter 
must  be  small  (see  ref.  1.1  cases  2  and  5). 


3*  Solution  of  the  Potential  Equation 


Assuming  harmonic  motion  of  angular  frequency,  v, 
equation  (2.2)  becomes, 

2 


where 


■  0 

"  a 


0  (x,y,z)  eivt  =  0(x,y,z,t) 


.(3.1) 

.(3.2) 


The  potential  0  is  subject  to  the  following 
boundary  conditions | 

(a)  0  is  bounded  everywhere  in  the  flow  and  at 

infinity  all  disturbances  should  disappear  in  the 
proper  manner,  thus,  _ 

i)  0>  0Z,  0~~>  0  as  \jy2  +  z 2 - 00 

ii)  The  solution  at  infinity  should  represent 
waves  travelling  outwards  from  the  origin, 

(b)  At  any  point  on  the  wing  the  prescribed  normal 
velocity  must  be  equal  to  the  normal  derivative 
of  0  at  that  point. 

Let  the  motion  of  the  point,  (x,y),-  on  the 
wing  be  represented  by, 

z  =  f(x,y,t) 

=  f  (x,y)  eivt  . (3.3) 

f  (x,y)  will  be  referred  to  as  the  *  deformation 
function.  * 


According  to  the  usual  assumptions  of  the 
linearised  theory  the  vertical  velocity  of  the  point 
(x,y)  is  given  by, 

dz  dz  TT  dz 

w  =  at  =  at  +  u£o  ax 

=  iv?(x,y)  +U;j0fx(x<y)  eivt  , 
or, 

* 

w(x,y)  =  iv  f(x,y)  +  U,0fx(x,y)  . ..(3.4) 

The  condition  is  satisfied  over  tie  projection  of  the 
wing  in  the  plane  z  =  0  and  takes  the  form, 

”'(x,y)  =  $z(x,y,0)  =  iv?(x,y)  +  ^fx(x,y)  ..(3.5) 

(c)  Outside  the  wing  and  outside  the  wake  0(x,y,z,t) 
must  be  continuous  in  planes,  x  =  constant,  and 
since  it  is  antisymmetric  in  z  must  satisfy  the 
condition, 

0(x,y,O,t)  =  0  , 


By  transforming  equation  (3*0  to  the  Elliptical 
Coordinates,  (£,ti) ,  where, 


y  =  s  cosh  n  cos  £  ") 

f 

z  -  s  sinh  ri  sin  £  ) 


(3.6) 


and  s(x)  is  the  local  semi  span,  Merbt  and  Landahl  (Ref  8) 
have  derived  a  solution  in  terns  of  Mathieu  functions. 


Using  che  notation  of  reference  18  the  solution 
takes  the  form,  ,f} 

£(x,ri,£)  =  'dL_.v  Pn  Ne^2^  (n,k)  sen  (g,k) 
n=1 


with  k  = 


.(3.7) 


*  The  use  of  a  'bar*  over  a  symbol  denotes  that  it  is  the 
amplitude  of  the  harmonically  oscillating  quantity  represented 
by  the  symbol  itself. 
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where  the  coefficients,  p  ,  are  to  he  determined  from 
"boundary  condition  (h),  n  equation  (3.5)  which  becomes 
in  elliptical  coordinates, 

^  (x#0,£)  =  s(x)  w  (x,s  cos  £)  sin  £  ..*...(3.8) 

Differentiating  (3«7)  with  respect  to  p  and 
putting  p  =  0  (on  the  wing )  gives. - 


^(x,0,S)  = 


writing. 


s'"  p 

— 7' 
n=1  • 

fW2> 

dp  v  n 

(n,k)j 

sen(^,k) 

p=0 

uo 

5a  Pn  sen(?,k) 

n=1 

• 

.(3.9)- 

• 

.(3.10) 

p=0 


Now  if  w(x,£)  is  bounded  and  is  a  continuous 
function  of  £,  the  series  representing  $  (x,0,£)  in  (3*9) 
will  be  uniformly  convergent,  ^ 

multiplying  (3.9)  "by  se  (£,k)  and  integrating  over 
the  range,  0  to  %,  the  coefficients  P  may  be  determined 
in  an  analogous  manner  to  Half -Range  PourSer  Coefficients 
since  an  orthogonality  relation  exists  for  the  Mathieu 
function  sen(£,k)  (see  ref.  16). 

The  Pn  are  thus  given  byj 


pn(k)  =  T"  |  sin  ^  sen(^,k)  <3^  (3.11 ) 


fi71 


U  o 


and  finally  the  from  (3.10), 

The  solution,  (3«7)»  is  now  completely  determined 
and  the  pressure  distribution  on  the  wing  will  be  given  from 
(3.7)  with  p  =  0, 

As  discussed  in  section  2  it  is  possible,  under 
certain  conditions,  to  suppress  the  time-  dependent  term  in 
equation  (3.1 )  and^tbe  solution  (3*7)  then  reduces  to,- 

$  =  e“nri  L  sin  n£ 

n=1  n 


(3.12) 
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where  the  L  are 
n 

given  by, 

-n\  = 

= 

fV7C 

“  j  w.sin 

sin  ng^  .  d£.j 

u  o 

m  i 

On  the  wing, 

J  --  • 

i  XI 

P’ 

_n 

n  * 

n=1 

sin  n£  . 

The  differential  pressure  across  the  vdng  plane 
is  given  as, 


Ap 


2p»V 


z=+0 


(3.15) 


4.  Symmetric  Flutter  Characteristics  Typical  of  a  Slender 
Delta  T/ing 

The  simplest,  pointed,  low  aspect  ratio  wing 
satisfying  the  assumptions  of  the  ’Slender  Body  Theory'  is 
the  slender 'delta  wing. 

Accordingly  the  analysis  as  developed  in  section  3 
is  applied  to  the  wing  shown  in  Figure  2. 

To  describe  the  possible  flutter  modes  of  the  wing 
a  reference  axis,  xQ(y)  is  used  (see  Fig.  2)  given  by  the 
equation, 

2c 

xQ  =  mcr  +  (1  -  \)  (1  -  m)  tyi  . (4.1 ) 

The  applicability  of  an  axis  such  as  this  to  delta 
wings  is  discussed  by  Y/oodcock  (ref.  17). 

For  any  particular  flutter  motion  it  is  then 
prescribed  that  sections  parallel  to  the  line  of  flight  will 
twist  about  the  reference  axis  according  to  some  modal 
function,  such  sections  remaining  themselves  undistorted, 
whilst  the  reference  axis  itself  translates  according  to 
another  modal  function,  each  degree  of  freedom  so  involved 
being  associated  with  a  Lagrangian  generalised  coordinate, 


+  Details  of  the  wing  are  given  in  Appendix  I'  . 


The  generalised  coordinates  are  defined  at'  reference 
sections  given  by 

|y!  =  1  . . (4.2) 

and  all  motions  are  measured  relative  to  the  mean  position  ; 
of  the  wing  (plane,  z  =  0). 

A  non-dimensional  spanwise  parameter,  6,  is 
introduced  such  that 

y  a  61  . . (4.3) 

and  }S|  =1  at  the  reference  sections. 

Each  degree  of  freedom  will  lead  to  an  equation  of 
motion  and  a  generalised  force,  Q  ,  which  can  be  expressed 
conveniently  in  terms  of  force  coefficients  as  (omitting  e^v"^) , 


* 


Q. 


P#> z 3u-2 


O'-' 


+  b 


rs 


iw  + 
m 


% 

..(4.4) 


on  the  assumptions  of  the  linearised  theory,  where, 

s  =  r  =  number  of  degrees  of  freedom 

and  <x>m  =  mean  frequency  parameter 

vcm  Note cm  is  the  mean  chord  of  the* 
*  half -wing.  See  Appendix  I. 

In  what  follows  the  suffices  +  and  -  will  indicate 
whether  a  function  applies  only  for  y  >  0  or  y  «<  0 
respectively. 


4.1 »  Uncoupled  Modes 

Let  there  be  one  uncoupled  mode  in  flexure  and  one 
uncoupled  mode  in  torsion  described  by  the  modal  functions, 
h(S)  and  H(6)  respectively,  so  defined  that, 

!h(+l)i  =  |h(+1)|  =  1 . (4.5) 

ivt 

If  ze  represents  the  translation  of  a  point 
on  the  reference  axis,  measured  from  the  mean  position, 


*  The  generalised  coordinates  and  forces  are  amplitude 
functions  but  the  bar  notation  is  not  used  in  their  case. 


-13- 


positive  upwards,  and  ae  the  rotation,  positive  leading 
edge  down,  of  the  section  through  that  point,  then, 


z  =  z 


,h(S)  -j 


and 


a  ~  aQ H(8)  Jj 


.(4.6) 


■where  z  and  a  are  the  amplitudes  at  the  reference 
o  o 

sections,  y  =  +  l. 

Now,  if  the  generalised  coordinates  are  chosen  so 


that, 


and 


z  \ 

...  _2  I 

l  [ 

°m  f 

^2  =  r  Go 


=  i 


.(4.7) 


then  the  deformation  function  f(x,y)  of  equation  (3.3)  takes 
the  form, 

...,.,......(4.8) 


I  -  Zhq.  +  (x-xq)  ~  Hq2 

m 


It  will  he  convenient,  for  the  aerodynamic  problem, 
to  consider  the  functions,  h  and  II,  each  to  be  polynomials 
in  6 1  thus,  for  symmetrical  modes} 


|  c  I  r  'j 
=  ^  |6| 


and 


H  =  S.G  Is) ' 

S  3 


.(4.9) 


r,s  ~  0,1,2,  ... 


Equation  (4.8)  becomes,  for  the  half -wing  for 
which  y  0, 

K? 


?+  =  1  j-^Sr6  (  + 


r 


—  =  ^>"  ct  5s  ( 


m  f  s 


. » . , (4.io) 


forces  and  Qg. 


It  is  now  required  to  find  the  two  generalised 
and 

Thus, 

I  i 

&P.  (6?+)  dy 


Q1  •  5<L,  = 


S. 


(4.11) 


*  'S’  indicates  that  the  area  of  integration  is  over  the 

wing  planform  for  which  0<  s(x)  <  y  only. 
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1  ’"“l 

j 

where  j  &  p  dx  dy  j  is  the  total  (incremental)  aerodynamic 
force  at  the  poini  (x,y)  on  the  wing  and  is  given  by 
equation  (3 .15)  when  the  velocity  potential  derived  for  the 
assumed  deformation  function  is  substituted. 

Prom  (4.10), 

-  1  s 

and  the  force,  Q^,  is  seen  to  be  built  up  from  a  sum  of 
integrals  of  the  form, 

Ap*  Sj.  .  dx  dS  . . ....(4.12) 

In  the  same  way  the  force,  Q^,  is  expressed  as  a 
sum  of  Integrals  of  the  form 

_  2  !>f' 

j  Q2|  a  “  j  (x-xM)  ^P,  Gb.  6s.  dx  d5  (4.13) 
-  -  s  m  JO 

It  will  be  clear  that  many  of  the  integrals  (4.12) 
and  (4.13)  will  be  identical,  apart  from  constant  factors. 


4.2.  Coupled  Modes 

The  deformation  function,  f ,  now  takes  the  form* 

f  =  «s=LUh(5)  +  (x-x  ) 

r  (,  m  J 

there  now  being  r  degrees  of  freedom. 

The  functions  h  and  H  are  as  defined  before  in 
equation  (4.9)  with  r  =  s  so  that  equation  (4.14)  becomes, 


?+  =  *>»  %.  j  »  +  <w<h>  a 


-a  l  sr 

:  r  /• 

m  J 


and  Q  is  given  by  I 

(T- 

Qr  •  =  Jj  (6?+)qr  ay  . 


...(4.15) 
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Finally,  from  (4.15) > 


«r-*‘  , 

1 


&  Sr  +  — 


L  _a\  e 


.  aj  8  .  ax  as 


.(4.16) 


This  is  a  sum  of  integrals  like  (4.12)  and  (4.13). 


5*  Calculation  of  the  Velocity  Potential  and  the  Generalised 
Forces  for  the  Assumed  Modes  of  Paragraph  4 

The  coefficients,  P  ,  in  the  series  representation 

of  $  are  given  by  (3.1 1)  and  since  f(x,y)  is  expressed 
in  polynomials  of  (8l  ,  and  hence  of  jyj  ,  integrals  of  the 
following  type  are  met  with! 


*  s 

cos  ^  sin  ^  sen(51,k)  d^  . 


.(5.1) 


Such  integrals  can  he  written  as  the  sum  of  integrals 
of  the  form 


j  sin  p^.  sen(^1,k)  <3^ 


-(5.2) 


Using  the  Fourier  Series  expansion  of  se  (£,k) 
integrals  such  as  (5.2)  become} 

C'J"  j  v 

(k)  j  sinr^,  sin  . . (5.3) 


\7hen  s  is  even  (equation  5.1 )»  the  limits  on 
(5.3)  are  0  to  %.  quite  straight  forwardly,  as  indicated 
"by  equation  (3.11  5  and  only  a  finite  number  of  terms  is 
obtained  for  (5.1 ).  Uhen  s  is  odd,  owing  to  the 
assumption  of  symmetry,  the  limits  on  (5.3)  reduce  to  0  to  %/2 
(or  %/2  to  %)  and  an  infinite  series  is  obtained  for  (5*3) » 
and  hence  for  (5.1  )•  However,  only  a  few  terms  need  be 
retained  in  practice. 

The  velocity  potential  on  the  wing,  is  now 

fully  determined  and  the  corresponding  loading  is  given  by 
equation  (3*15). 

The  generalised  forces  give  rise  to  integrals  like 
(omitting  constants), 


6r  dx  d6 


.(5.4) 
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These  integrations  must,  in  general,  "be  done  by  a 
graphical  or  numerical  means,  except  when  the  time-dependent 
term  in  the  potential  equation  (3*6)  is  suppressed  and  the 
Mathieu  functions  take  on  their  degenerate  forms. 

The  application  of  the  analysis  to  antisymmetric 
flutter  modes  follows  the  same  general  lines  as  given  for 
symmetric  modes. 


5.1 .  Equivalent  Constant  Flutter  Derivatives 

By  analogy  with  the  flutter  derivatives  of  two- 
dimensional  (strip)  theory  it  is  possible  to  define  a  set 
of  'equivalent  constant  derivatives'. 

These  derivatives  are  constant  over  the  span  of 
the  wing  and  give  the  correct  generalised  forces  when  inter¬ 
preted  in  the  conventional  sense. 


The  lift  and  moment  on  a  strip  of  unit  width  are 
defined  in  terms  of  derivatives  such  as 


(ref.  is) 


where  the  'stiffness'  derivatives  include  the  'inertia' 
derivatives,  lg  ,  m^,  ,  ,  m^.  . 

Equivalent  constant  derivatives. 


(l  ) 

'  z* 

,  dj) 

i 

rs 

rs 

etc. 

(m  ) 

,  (m2) 

» 

rs 

rs 

are  defined  from 
Appendix  II. 

the  force 

coefficients 

As  with  the  force  coefficients  the  first  suffix 
refers  to  the  generalised  force  and  the  second  to  the  mode. 


*  See  Appendix  II  for  discussion  of  sign  convention, 
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Apart  from  the  analogy  with  '  two-dimensional  * 
derivatives  the  concept  of  equivalent  constant  derivatives 
is  useful  in  that  it  facilitates  direct  comparison  of  sets 
of  derivatives  derived  for  different  modes.  For  example, 
rect  derivatives  in  one  freedom  are  mode  independent  of 
modes  in  other  freedoms.  (See  Appendix  III). 


6.  Representative  Results 

The  preceding  analysis  has  been  applied  to  the 
case  of  a  triangular  wing  (Fig.  2,  k  =  0)  using  uncoupled 
modes. 


Equivalent  Constant  Derivatives  have  been  calculated 
for  the  flexural  modes} 


h(s)  =  !&lr 

and  torsional  modes} 

H(8)  =  !  6  f s 


r  =  0,  1,  2. 


s  =  0,  1. 


Modes  such  as  these  have  been  taken  in  pairs,  one 


in  flexure,  |5|r1  and  one  in  torsion, 
’sets’  of  derivatives. 


iSl 


si 


giving  six 


The  accompanying  table  of  numerical  results  (Table  i) 
shows  the  order  of  the  derivatives  and  their  signs  (for  m  =  -g-) 
and  a  set  of  general  expressions  for  the  derivatives  is  given 
in  Appendix  III  together  with  the  results  of  a  calculation  on 
a  cropped  delta  for  r^  =  s^  =  0  only. 

The  'damping'  derivatives  ,  m^  and  m^  are 
plotted  against  'm'  in  Figs.  8  $  and  10 

By  taking  r^  =  s^  =  0  and  m  =  1  the  derivatives 

are  obtained  for  a  rigid  pitching  and  plunging  v/ing  referred 
to  the  trailing  edge  -  use  of  the  usual  transformation 
formulae  then  refers  the  derivatives  to  any  other  axis. 

This  has  been  done  for  a  triangular  wing  (X  =  0), 
a  cropped  delta  (X  =  1/7)  and  a  rectangular  wing  (X  =  l) 
for  an  axis  at  0.500  c  and  the  results  are  presented  in 
Table  II.  r 

In  Fig.  11  the ’cross-damping’  derivatives  m^  , 
have  been  plotted  against  ’m'  for  these  wings. 

In  this  case  of  a  pitching  and  plunging  wing  the 
generalised  forces  and  Q2  have  simple  interpretations 


and  do  in  fact  represent  the  unsteady  lift  and  moment 
amplitudes  on  the  complete  wing,  i.e,. 


p  U?..  S  p  U2  Z-5 


x  -  ,  (+ve  upwards) 

cr 


M  Q2  .  -iZ2  /  „  ■  ,  n 

—5 —  =  — — 5 — =•  x  — s*  ,  (+ve  nose  dounvards) 


2  “  2  ^ 
p  U7  Sc  p  U  „  V 

rCv  (■*  r  •v-,  ,>5 


,(6.1) 


The  expressions  for  lift  and  moment  vd.ll  be  in 
terms  of  the  dimensionless  amplitudes. 


f  Z  V 

(t) 


and  a 


and  the  relevant  frequency  parameter  will  be. 


v  r  /  2  \ 

U,/;j  =  °m  Ij+X/  ‘ 


It  is  convenient  in  this  connection  to  define  a 
set  of  force  coefficients  for  rigid  motions  only  since  in 
later  paragraphs  unsteady  lifts  and  moments  on  rigid  bodies 
and  wing-body  combinations  are  considered. 


Coefficients  Lz>  L,  ...  etc.  are  defined  by  the 
expressionsj 

T  /  z  \ 

- =—  =  -  (L  +  iw  L.)  I  —  j  -  (L  +  ica  L.)  a 

.  tt2  o  z  r  z'  \  o'  'a  r  fl'  o 

and  ..  •  z  N 

“V,  =  »az  +  ^ j  +-  <i!a  *  V  “o 

^  . ....(6.2) 


and  these  will  be  referred  to  as  'Rigid*  force  coefficients. 


As  for  the  definitions  giving  the  equivalent 
constant  derivatives  (Appendix  II)  these  rigid  farce  coeff¬ 
icients  are  signed  to  agree  with  the  normal  British  flutter 
sign  convention 
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TABLE  I 


-  K  w1 

+  la  (A)"1 

n,  r 

s  X 

0 

■ 

1 

2 

HB 

0 

1 

2 

0 

1 

.0296 

.0296 

.00768 

.OO768 

.00374 

.00374 

0 

1 

.781 

1,00 

.996 

.955 

ma  ort1 

■s 

0 

1 

2 

r 

s  X 

0 

1 

2 

0 

1 

.00431 

.00214 

,00450 

.00200 

0 

1 

.0530 

.0334 

.0530 

.0334 

.0530 

.0334 

+  l&  (A)-1 

+  ia  OR)"1 

9W8P 

0 

1 

2 

SS 

0 

1 

2 

0 

i 

.786 

.786 

.954 

.954 

.981 

.981 

0 

1 

.735 

.915 

1.02 

.905 

1.30 

.817 

m5  (AR)"1 

-  (51  r1 

Si 

1 

2 

ss 

0 

1 

2 

0 

1 

.0542 

.0900 

-.0223 

.0379 

*.135 

-.0405 

0 

1 

.114 

.195 

_ 

.114 

.195 

.114 

.195 

Derivatives  for  m  =  ■§■;  =  0,2 


Equivalent  Constant  Derivatives  for  rigid  wings  plunging  and 
pitclring  about  an  axis  at  -jr  c 


Equivalent 

Constant 

Derivatives 

Triangular 

wing 

(1=0) 

Cropped 

Delta 

(1=1/7) 

Rectangular  i 
Wing 
(1=1) 

l 

a 

at  2 

- 

(.524) 

-.585  »m2 

0t  2. 

"  4um 
(.785) 

l5 

+  4 
(.785) 

+.785 

at 

+  4 
(»785) 

l 

a 

%  2  % 

-  T2  wm  *  4 

(.262)  (.785) 

-.2450^  +*785 

*  f 

(.785) 

l. 

a 

♦  ft 
(.961) 

+1  *07 

+ 

+  8 
(1.18) 

m 

a 

at  2 

12  wm 

(.262) 

2 

.245  o>m 

0 

at 

(.197) 

-.0845 

7C 

8 

(.392) 

m 

a 

Ot  2  0t 

20  Wm  "  T£ 
(.157)  (.197) 

.13?(4  -.0845 

7t  2  ,  0t 

48  “m  +  8 
(.0655) (.392) 

at 

'  8 

j  (.392) 

-.386 

at 

“  15 
(.197) 

NOIEJ 


Figures  in  brackets  are  decimal  equivalents  of 
fractions  of  at. 
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lender 


of  Revolution 


A  set  of  flutter  force  coefficients  for  a  pointed 
slender  body  of  revolution  can  be  calculated  in  an  analogous 
manner  to  those  of  the  wing  by  adopting  polar  coordinates 
instead  of  elliptical  coordinates  when  solving  the  potential 
equation  and  in  the  specification  of  the  boundary  conditions. 

Probably  the  only  case  of  interest  is  the  rigid 
pitching  and  plunging  body  and  accordingly  this  case  will 
be  dealt  with.  The  cartesian  coordinate  system  for  the 
body  is  the  same  as  for  the  wing  and  is  shown  in  Figure  3 

In  each  cross-section,  x  =  const.,  tale  polar 
coordinates! 


y  =  r  cos  £  ^ 
z  =  r  sin  £  ( 


.(7.1) 


then  the  potential  equation  (3*6)  transforms  toj 


1  <L 

r  *  ar 


(rf)  *' 


i-  +  l 

r2  as2  a2 


0 


.(7.2) 


Consider  the  body  movements  to  consist  of  a 
vertical  translation  (+ve  upwards)  and  pitching  about  an 
axis. 


X  =  XQ  = 


“b  lB 


.(7.3) 


parallel  to  the  y-axis  (nose-down  pitching  +ve). 

By  analogy  with  equation  (4.6)  we  define  zq  to 
be  the  amplitude  of  the  displacement  of  the  point, 
x  =  xq,  on  the  body  axis  and  aQ  to  be  the  amplitude  of 

the  inclination  of  the  body  axis  to  the  Ox  axis.  Then 
the  motions  produce  at  a  point,  x,  on  the  body  axis  the  total 
upward  displacement, 

(z  +  (x-x)  a )  eivt  . 

Taking  the  body  length,  lB>  as  a  reference 
length  the  vertical  velocity  (equation  3*4)  is, 

”  iv  1  hi  iz)  +i  -o  +  Z)  °J  . (7,4) 

i  '•  '•  .} 

Writing  V  (x)  for  the  local  cross-sectional  area 
of  the  body,  the  potential  near  the  body  takes  the  foxraj 
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(&)  =  "  Tf(x>  .  sin  ^  .f . .(7.5) 

.  r  -*  R  ^ 


"body  is 


where 


The  unsteady  lift  and  moment,  L  and  M  follow 
by  integration  of  (7.6)  along  the  body. 

The  'rigid*  force  coefficients  of  equations  (6,2) 
have  been  calculated  for  a  cylindrical  body,  with  a  conical 
nose  as  shown  in  Fig.  4.  These  are  based  on  the  Aspect 
Ratio  of  the  geometrically  similar  wing  having  its  root-chord 
equal  to  the  length  of  the  body  and 

b  R 

2  =  ~  (see  Figure  7) 

where  Rq  is  the  maximum  (base)  radius  of  the  body. 

The  coefficients  are  given  in  Appendix  IV. 

By  putting  cr  =  1  and  t\  =  (1-X)  it  will  be 
seen  that  the  expressions  of  Appendix  IV  are  identical  with 
those  that  would  be  given  fear  the  rigid  cropped  delta  wing 
using  the  equivalent  constant  derivatives  of  Appendix  III 
with  B  =  0  and  equations  (6.1 ), 

8.  The  Slender  Ting-Body  Combination 

8,1  The  rigid  pitching  and  plunging  combination 

A  set  of  'rigid1  force  coefficients  will  now  be 
derived  for  the  slender  wing-body  combination  shewn  in 
Figure  8, 


The  pressure  at  ary  point  on  the  surface  of  the 


-  TT  (  3 II  xv  „ )  sxn  £ 

P  =  p  ..  U  o  1  x —  +  tt  .II >  - a 

h.W  v-'-'  i  d  x  U,?v  {  to: 

Il(x)  =  -  w(x)  T(x) 


.(7.6) 


This  problem  will  be  dealt  with  rather  differently 
from  the  wing  and  slender  body  cases  in  that  the  velocity 
potential  will  be  found,  not  directly  as  a  solution  of 
Laplace  *  s  Equation,  but  from  the  two-dimensional  potential 
for  incompressible  flow  normal  to  a  flat  plate. 

The  required  potential  will  not  generally  satisfy 
the  two-dimensional  wave  equation  (2,2)  and  hence  the 
solution  will  be  subject  to  similar  restrictions  as  the  wing 
solution  for  k  *-••->  0. 
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Using  the  Joukowski  Transformation  (see  Fig,  6) 
the  velocity  potential  for  the  flow  around  the  body  config¬ 
uration  of  Fig,  6  (see  ref,  4)  due  to  a  motion  w  of  any 
section  in  a  fluid  at  rest  oan  easily  be  found  and  conven¬ 
iently  expressed  in  two  parts, 


where  B(0)  is  the  potential  on  the  body  (r=R)  and  w(0) 
is  the  potential  on  the  wing  (£  =  0  or  x,  y  =  r), 


It  wall  be  clear  that  the  force  coefficients  for 
the  wing-body  combination  of  Fig.  5  can  be  considered  to  be 
the  addition  of  two  sets  of  force  coefficients;  viz., 

(i)  the  force  coefficients  for  a  triangular  wing  on 
a  cylindrical  body,  downstream  of  the  lateral 
plane  through  the  wing  leading  edge  and  body  junction. 


(ii)  the  force  coefficients  for  a  pointed  body  upstream 
of  the  wing  leading  edge. 

The  coefficients  (ii)  have  been  calculated  in 
section  7  (Appendix  IV). 

The  coefficients  (i)  can  be  calculated  using  (8.1 ) 
with  the  axes  and  notation  of  Figure  7* 


For  the  combination,  the  velocity,  v,  will  take 
the  same  form  as  for  the  body  alone,  i.e.  equation  (7.4), 
thus. 


5(x)  =  ivi  O  If]  -  Vo  +  W“o 

t> _  \  r’ 


+  «0*j  (8*2) 


The  loading  distribution  is  given  by  the  pressure 
equation  (3 .15)  and  lift  and  moment  byj 


-;or  f\* 


L  = 


1  i  (A\p)B.8y 


J 


M  =  ^ 


r-°r  f  |'E 


-Ip 


i 

(&p)w.ayj  8x 


O  otc 


jUo 


(Ap)B.8y  + 


R 


3  “'n 

i 

(Ap)w.<3y|'  (x-xo)  dx 

R  ......  i . (8.3) 


and 
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Complete  expressions  for  lift  and  moment  on  the 
clyindrical  body  and  triangular  vring  lead  to  the  ’rigid1 
force  coefficients  which  are  given  in  Appendix  V,  These 
coefficients  like  those  of  the  body  are  based  on  the 
triangular  wing  having  root  chord,  c  ,  and  maximum 
semispan  sQ  =  b/2  =  R  r 

The  force  coefficients  for  the  whole  wing-body 
combination  of  Figure  5  are  given  in  Appendix  VI  and  the 
variation  of  the  ’damping*  force  coefficients  L,  ,  L,  ,  M.  , 
Mg  with  a4  is  shown  graphically  in  Fig.  12.  ' 

In  adding  the  appropriate  coefficients  of 
Appendices  TV  and  V  the  definitions  of  Figure  5  and  Appendix 
VI  were  used  and  again  the  triangular  wing  is  used  as  a  basis. 


9.  Discussion 


The  use  of  the  ’Slender  Body  Theory’  for  unsteady 
flow  problems  leads  to  a  solution  for  the  aerodynamic  forces 
which  dees  not  involve  long  computation  and  many  geometrical; 
and  other  parameters  can  be  carried  along  in  the  analysis 
without  having  to  be  specified  definitely  at  the  outset. 

The  restrictions  of  the  theory  as  discussed  in 
section  2  seem  to  be  somewhat  severe  but  there  is  evidence 
to  show  (ref,  8)  that  far  a  rigid  triangular  wing  of  aspect 
ratio,  1 ,  at  a  Mach  number  of  1,25  and  far  a  frequency 
parameter,  to  ,  up  to  6,  the  theory  appears  to  be  quite 
valid.  Furthermore,  results  far  an  aspect  ratio  of  0.5  show 
that  when  the  time  derivative  terms  are  neglected  the  results 
differ  from  those  given  by  the  complete  solution  only  if 
w  >  2  for  a  Mach  number  range  of  0  -  1.25  (rigid  triangular 
wing). 


Owing  to  the  reed  to  evaluate  several  terms  of 
the  Mathieu  function  series  when  deriving  the  full  solution 
it  is  much  longer  than  the  simplified  case  (for  snail  root 
frequency  parameter)  and  it  would  always  be  worthwhile  to 
question  whether  the  full  solution  is  really  necessary  in 
any  specif  led  case. 

With  the  type  of  wing  to  which  this  analysis  can 
be  applied,  it  is  very  unlikely  that  the  root  frequency 
parameter  will  exceed  about  0,5  so  that  in  many  cases  the 
simplified  approach  would  suffice. 

The  force  coefficients  Y,h,c  of  equation  (4.4) 
are  dependent  on  Mach  number  and  frequency  only  through  the 
parameter,  k,  in  the  general  solution,  consequently,  in 
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the  simplified  case  which  implies  that  k — 0,  the  coeff¬ 
icients  are  independent  of  frequency  and  llach  number. 

It  has  been  found,  both  experimentally  and 
theoretically  that  the  variation  of  flutter  force  coefficients 
with  frequency  decreases  as  aspect  ratio  decreases  so  that 
this  is  not  a  surprising  result  from  a  theory  which  is  correct 
far  M — &  0. 

The  preceding  remarks  can  be  taken  to  apply  equally 
well,  in  principle,  to  the  wing-body  combination. 

It  is  interesting  to  note  that  the  analysis  used 
by  Lawrence  and  Gerber  (ref.  19)  (subsonic)  when  taken  to 
the  limit  M  — 0  gives  results  for  a  rigid  wing  which 
agree  with  those  found  here  and  by  Garrick  (ref.  1 4) . 

In  this  connection  it  is  also  interesting  to  study 
their  results  when  plotted  against  aspect  ratio.  The  slopes 
of  the  curves  (force  coefficients)  at  zero  aspect  ratio  are 
correctly  those  given  by  ’Slender  Body  Theory'  but,  in 
general,  the  curves  depart  from  their  original  tangents 
extremely  rapidly.  It  might  be  suggested  therefore  that 
force  coefficients  derived  using  'Slender  Body  Theory',  if 
applied  outside  their  range  of  reasonable  validity,  will 
give  magnitudes  which,  in  general,  will  be  very  different 
from  the  'true'  values.  (See  also  refs.  2o  and  21  for  which 

M  =  3). 


Fig,  10  shows  that  for  an  axis  at  the  trailing 
edge  of  the  wing,  no  matter  which  torsional  modes  are  chosen, 
the  direct  damping  derivative  m^  is  zero  indicating  that 
an  undamped  pitching  oscillation  would  be  possible  -  for  all 
axis  positions  0  <.  m  <  1  the  derivative  gives  positive 
damping, 

Ythen  w  —  the  'rigid'  force  coefficients  give 

the  values  of  lif?  and  moment  for  the  steady  case  as  found 
by  Jones,  Spreiter  and  others  (refs,  7  and  4.) 
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AHENDIX  I 


The  Cropped  Delta  - 

Definitions  and  Geometrical  Properties 

See  Figure  2  ;  - 

Mean  chord,  cm  =  ■§•  (l  +  X) 


-Area,  S  =  h  c 

1  m 

h2  h 

Aspect  Ratio,  M  ~  “  c” 

m 


(D 

(2) 

(3) 


local  semispan,  s  = 


.x  for. 


2c  (l-X) 

r  OCx<(l-X)c 


U 

Reference  section  position,  l  =  -g 
Spanwise  parameter,  8  =  ^ 


2c 


local  chord,  c  =  - -  •  (1  -  t(1-X)&) 


(1+X) 


Ratio, 


l  JR  ,x 


m 


2c 


Reference  Axis,  xq  =  mc^  +  (1”X)(l**m)  jyi 

=  A  +  B  \b\ 


(4) 

(5) 

(6) 

(7) 

(8) 

(9) 


thus, 

and 


A  _  2m 
cm  1+X 

ss  .  k  h=hl  (1.m) 

Iom  (1+X) 

VC 

,  ra 

Frequency  parameter,  to  = 

(Mean)  V' 


(10) 


(11) 
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AHENDIX  U 


Equivalent  Constant  Derivatives  »  Definitions 


The  following  equivalent  constant  derivatives  are 
appropriate  to  a  system  with  the  usual  British  sign  convention 
i.e.  z-axis  downward,  lift  positive  upwards,  moment  and  angle 
of  attack  positive  nose-up* 

(a)  Uncoupled  Modes 

i.  One  mode  in  flexure ,  one  mode  in  torsion 

(*iA 


+  c. , )  =  l 

'  l11  n  11 '  v 


^12<“b  +  012>  =  \ 


H>1°£  +  °21 )  =  <■%> 


h2  .  as 


viA 

o 

•jiA 


h.H  d8 


m 


VO 

yiA 


—  H.h  d5 
c 
m 


OV-t  2 

(-*22%  +  °22>  =  ]  (^)  36 


■vi  A 


b  w  =  0)  .  1. 

11  m  m2 


b12  <%  =  “m  !a 


‘'  f  O 

•vi  A 


—  .  hfc  d8 
c 
m 


tio 


(t) 


h.H.  d8 


b21  "b  =  “b  (-ml) 


\lA  2 

H.h  d6 

0 

ra/ 


b22  “m  =  “n  (-“a) 


U  o 

«lA,.  -  3  , 

\  IT  d6  . 


1  (k) 

V  O  v  ' 


ii,  r  modes  in  flexure*  s  modes  in  torsion 
Derivatives  such  as  ( lz)  ^  and  (mz)  ^ 


rr 


sr 
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Yri.ll  be  defined  in  the  same  manner  as  in  (i)  by  integrals 
such  as, 

fliA  ■  f)i  A 

l  h  .  h  ,  d6  and  —  H  .  h  ..  d6  respectively. 

J  o  U  o  m 


(b)  Coupled  Modes 


(-Y  w2  +  c)  =  (I  )  h,hd5 

N  ‘rs  m  rs,  N  z' _  r  s 


rs  ,  i 


(~Y  co+c)  =  ( l  )  —  h  H  d5 

v  ‘rs  m  rs/.  '  aJ  c  r  s 

I  rs  ,  |  m 


h  H  d5 


rs  I  v  m 
•J  o 


!  c  I  r  s 


etc.  -  by  analogy  with  (a)i, 

2 

A  term  such  as  (-y  w  +  c  )  is  obtained  from 
v  ‘rs  m  rs  j 

z 

the  real  part  of  the  term  in  the  total  expression  for 

which  involves  both  h  and  h  j  and  (b  w  )■,  .  ,  ,  .  - 

r  s'  ^  rs  rlj  is  obtained 

from  the  imaginary  part  of  the  sane  term:  other  terms  are 
obtained  in  a  similar  manner. 
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(a)  Equivalent  Constant  : 


3  (.516 


-  ,349m)  R 
3  (.139  -  .205m)  w2  .  R 
3  (.182  -  ,349m) 

|^I(.142  -  .¥Xta  +  ,317m2)  w2 

(  ) 

+  (-  .121  +  .288m)) 

J 

§  (-  .353  +  .703m  -  ,354m2) 

1?8  *  wm  m 


f? 

•  & 

6%  < 
1 

\k  <-430 

k 

-  ,690m) 

0 1  -  .0625,1  M 

”  J 

15tc 

ir 

■  (.690  - 

•49Cta)  JR 

3ft 

2 

(.082  -  . 

115m)  a)2 
m 

R 

IJft 

4 

■  (.0730  - 

,123m)  R 

321  1 
4  1 

i  2( ,0711  ■ 

-  ,200m  + 

•158m2) 

2 

°ra 

‘* 

-  (.121  - 

,288m)  .1 

R 

%  (.178  -  ,354m  +  ,176m2)  R 


3  { ■§■(•220  -  .353m)  u>2  “  -333  >  £ 

L  J 

3  (.2*42  -  .2724m)  JR 

3  (.HO  -  ,176m)  w2  ,  M 

3  (.107  -  .271tfn)  ^ 

1 8  (.00359  -  ,0100m  +  ,00800m2)  w2 
‘  -  ^  (.0430  -  .0930m )l  M 

30  (-.0250  +  ,02f92m  -  .0243m2)  S. 


3  2  _ 

r-  .  co  iR 
m 


(.0592  -  .0911m)  co2  -  1 1 
|  (0.540  -  0,355m)  M 


6  (.0395  -  .061 2m)  </  .  & 


15-  (.052(0  -  ,118m) 


15.  |  (.0287  -  .0800m  +  ,064Cto2)  co2 

*“  -  ^  (.0430  -  ,0930m).!  R 

%  ) 

30  (-  .0250  +  ,0492m  -  .0243m2)  R 
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r  =  2,  s  =  1 . 


5%  2  _ 

TZo  W  ,  M 
168  m 


-  jf  (.0792  -  .117m)  to2  -  l.ooj 
+  (.206  -  ,194m)  JR 


& 


ra  = 


m. 


m 


m. 


-  40  (-  .004-95  +  .00735m)  (/  JR 


30  (-  .0359  +  .0691m)  JR 

£  <  (.0287  -  .0800m  +  .06243m2)  w2 
2  2  -  m 


-  i  (.0430  -  ,0930.1 


jr 


30  (-  .0250  +  .0492m  -  ,0243m  )  JR 


(b)  Equivalent  Constant  Derivatives  for  Cropped  Delta 
(X,  general) 
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b  =  lr.±^  -  3/2->2i 

9  15(1  +  X)2 

o  1  +  4X 
P10  =  20(l+X) 


Expressions  far  the  above  derivatives  when  X  =  l/7 


=  -  .585  </  .  JR 


m 


♦  5  * 


1 ( •  i  4-1  -  .226m)  o) 


m 


-  I 

“H 


JR 


lft  =  +  .2i  1 7c(i  .86  -  1.12m)  JR 

mz  =  7t(.141  -  »226m)  w2  JR 

=  -  .845'*  (.0944  -  .281m)  JR 

=  ,845k  I  (.159  -  .384*  +  .284m2)  co2 

a  ^ 

-  (.0944  -  .28 moj  iB 
=  ,645k  (-  .325  +  .650m  -  .325m2)  JR 


are 


•  • 
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AFEENDIX  IV 


'Ric id’  Force  Coefficients  for  a 


with  conical  nose 


(see  Figure  4) 


Definitions I 


Length  of  "body  =  l B 

Length  of  conical  nose  =  A  i-g 

Reference  axis  at  distance  nig  from  nose 

Frequency  Parameter  oig  =  v  , 

The  force  coefficients  are  based  on  the  wing  having 
a  geometrically  similar  planform  to  the  body,  thus; 


Rn  =  9  •  and  Z-j  =  c  , 


so  that 


(«■(*) 


and  on,  = 
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APEENDIX  1/T 


*  Rigid1  Force  Coefficients  for  a  Slender  V/ing-Body  .Combination 

(see  Digure  $) 


Definitions,  -  ’ 

Total  length  of  combination  =  Vc 


l 


(1) 


B 


Ratio  of  body  length  to  wing  root  chord  =  —  =  T 

c 
r 


(2) 


But  Vc^  =  ( ig  +  0^(1 -cr)  )  ,  see  figure  5, 
so  that  frcm  (2)P 

T  =  V  -  (1  -  o')  .  (3) 


It  is  essential  that  the  position  of  the  reference 
axis  should  be  unique  when  measured  from  the  apex  of  the  wing 
(mc  cr)  and  from  the  nose  of  the  body  (nig  ,  ig).  This 

requires  that  (see  figure  5)  ! 


“B  lB 

=  »b  -  %  +  ”0°r  • 

(o'-m) 

i.e. 

mg  = 

1  -  c 

1  T 

(6) 

V  -  1+m 

or 

*B  = 

c 

V  -  1  +  cr 

(7) 

Length  of  conical  nose  =  /\  ig  =  Alc^  (8) 

This  gives  nose-length  as  a  constant  proportion  of  body  length 
ahead  of  wing  root.  If  a  nose  length  which  is  a  constant 
proportion  of  total  combination  length  is  stipulated  then  /  > 
must  be  replaced  by  an  expression  of  the  form; 


=A(  y  j  =  const. 

(9) 

Note! 

c  3 

In  the  above  definitions!  T  4 

0 

Y  4 

( 1— cr) 

(10) 
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